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We investigate remote surface scattering (RSR) by the Si02/Hf02 interface in Fully-Depleted Silicon-on- 
Insulator (FDSOI) devices using Non-Equilibrium Green’s Functions. We show that the RSR mobility is 
controlled by cross-correlations between the surface roughness profiles at the Si/Si02 and Si02/Hf02 inter¬ 
faces. Therefore, surface roughness and remote surface roughness can not be modeled as two independent 
mechanisms. RSR tends to enhance the total mobility when the Si/Si02 interface and Si02 thickness profiles 
are correlated, and to decrease the total mobility when they are anti-correlated. We discuss the implications 
for the high-ft/Metal gate technologies. 


High-K materials such as Hf02 have been introduced 
in Metal-Oxide-Semiconductor Field Effect Transistors in 
order to keep a tight electrostatic control over short chan¬ 
nels while limiting gate leakage currentsii They usually 
remain separated from the channel by a thin interfacial 
layer (IL) of Si02. Yet the introduction of high-K mate¬ 
rials leads to a systematic decrease of carrier mobilities, 
which is very significant at weak inversion densities, but 
can persist in the strong inversion regime 

Different scattering mechanisms have been put forward 
to explain this degradation. There is no doubt that 
charges trapped at the Si02/Hf02 interface or in the 
Hf02 layer (Remote Coulomb scattering (RCS)) make 
a major contribution at weak inversionj^^— Remote scat¬ 
tering by polar optical phonons (RPH) in Hf02 is also a 
serious candidate, especially in thin IL devices.— Much 
less attention has been given up to now to scattering by 
roughness at the Si02/Hf02 interface or equivalently by 
IL thickness fluctuations. At variance with RCS, this 
mechanism, known as remote surface roughness (RSR) 
scattering, is expected to be dominant in the strong in¬ 
version regime. 

RSR has first been investigated with semi-classical 
Kubo-Greenwood approaches in a different context, 
namely roughness at the Si02/Gate interface in polysil¬ 
icon gate technologies— The models were later ex¬ 
tended to Hf02 and Hf02/IL gate stacks d^d^ However, 
the different interfaces were assumed uncorrelated, and 
surface roughness (SR)/remote surface roughness mod¬ 
eled as two independent mechanisms. 

In this letter, we use Non-Equilibrium Green’s Eunc- 
tions (NEGF) methodsi^dd to investigate RSR in the 
latest Fully-Depleted Silicon-on-Insulator (FDSOI) thin- 
film technologies with high-At/Metal gates. We show 
that RSR scattering is dominated by roughness at the 
Si02/Hf02 interface, indeed prevails at large carrier den¬ 
sities (at variance with RCS), and decreases exponen¬ 
tially with IL thickness. More importantly, we demon- 
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strate that the RSR correction is controlled by cross¬ 
correlations between the Si/Si02 and Si02/Hf02 inter¬ 
faces. SR and RSR can not, therefore, be modeled as 
two independent mechanisms. RSR is beneficial when 
the Si/Si02 interface and the IL thickness profiles are 
correlated, and detrimental if they are in anti-correlated. 

Let us first discuss the basics of remote surface scatter¬ 
ing. The main effect of RSR is to modulate oxide thick¬ 
nesses hence the effective field within the semiconductor. 
This gives rise to potential and charge fluctuations re¬ 
sponsible for extra carrier scattering— “-12 In this picture, 
we expect stronger scattering from the Si02/Hf02 inter¬ 
face than from the Hf02/Metal gate interface, since the 
latter is usually much farther from the channel and is 
efficiently screened by the Hf02 layer. This will indeed 
be confirmed by our numerical simulations on FDSOI 
devices. From now on, we hence focus on RSR at the 
Si02/Hf02 interface. 

We can draw the main trends and identify the rel¬ 
evant structural parameters from semi-classical pertur¬ 
bation theory. Let A(r) and A'(r) be the SR profiles 
at the Si/Si02 and Si02/Hf02 interfaces, and i5(r) = 
A'(r) — A(r) be the variations of the thickness of the 
IL. According to Fermi Golden Rule, the scattering rate 
between states |k) and |k') with in-plane wave vectors 

k and k' is proportional to |(k'|iJF>|k)|^, where is 
the Hamiltonian of the SR-FRSR disorder, the over-line 
denotes an average over configurations, and we have dis¬ 
carded band indexes for simplicity. This matrix element 
can be expanded in powers of A and S: 

(k'|7JR|k) = A(q)A(q) -F B{q)5{q) -f ..., (1) 

where q = k' — k, A(q) and 5(q) are the Fourier trans¬ 
forms of A(r) and ^(r), and A(q) and B{q) are complex 
numbers. Hence, to second order in A and 6, 

l(k'|i?R|k)|" = |A(q)nA(q)|^ + |H(q)p^^ 

-F A*(q)H(q)A*(q)(5(q) -Fc.c. (2) 

The first term will be referred to as “pure SR” scat¬ 
tering (SR scattering is, indeed, usually compute d^^d^ 
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with assumptions best compatibleii with 5 = 0). It is 
proportional to |A(q)|^, which is the Fourier transform 
of the auto-covariance function Fa(R) = A(r)A(r -|- R) 
of the SR prohle at the Si/Si02 interface (r and R 
being in-plane vectors). The second term is “uncorre¬ 
lated” RSR scattering by fluctuations of the IL thick¬ 
ness. It is, likewise, proportional to the Fourier transform 
of the auto-covariance function ^^(R) = 5(r)5(r + R) 
of these fluctuations.— ““ The third (and fourth complex 
conjugate) terms are “cross-correlated” contributions to 
RSR scattering. They are proportional to A*(q)5(q), 
which is the Fourier transform of the correlation func¬ 
tion G'(R) = A(r)5(r -1- R) between the SR profile at the 
Si/Si02 interface and the IL thickness fluctuations. They 
have not been accounted for in Refs. [Ml- 

In general, one expects |R(q)| <C |^(q)|, so that the 
cross-correlations prevail over the uncorrelated term un¬ 
less G'(R) 0 - which is unlikely in the thin ILs where 

RSR might be signihcant. The above analysis therefore 
shows that SR and RSR can not be modeled as two inde¬ 
pendent mechanisms, as assumed in previous literature. 

We consider an exponential (auto-)correlation model 
for the Si/Si02 interface profile A(r) and IL thickness 
variations 5(r)— ^ 

Fa(R) = A2/(R) ; Fs{-R) = S^fiR) 

G{R.) = ASf{R)g, (3) 

where f{R) = and = 1.3 nm everywhere. 

g characterizes the cross-correlations (or similarity) be¬ 
tween A(r) and 5(r) (|(;| < 1). This is equivalent to 
the following model for the Si 02 /Hf 02 interface profile 
A'(r): 

Fa'(R) = A'(r)A'(r-f R) = A'V(R) 
iL(R) = A(r)A'(r + R) = AA'f{R)h , (4) 

where A'^ = A^ -|- 5^ -|- 2A6g and h = (A + 6g)/A' is 
the overlap between the two interfaces. They are paral¬ 
lel if h = 1, uncorrelated if h = 0, and “anti-parallel” 
if h = —1 [Conversely, 6^ = A^ -|- A'^ — 2AA'h and 
g = {A'h — A)/S\. Although the parameters A' and h are 
more intuitive;^ the equivalent (5, g) set is more prac¬ 
tical. Indeed, pure SR scattering simply corresponds to 
5 = 0; Also, this set provides a straightforward criterion 
for the integrity of the Si02 layer: The surface covered 
by holes (5(r) < —fm) in an IL with nominal thickness 
tiL is < 2.3% as long as 5 < tiL/2— ^ The present results 
and conclusions might not, therefore, hold if 5 tiL/2. 
The squared matrix element Eq. © thus reads: 


|(k'|RR|k)|^=/(q){|A(q)pA2 

+ 25Re[A(q)R*(q)]A5g} , (5) 

where we have neglected the oc |R(q)p5^ term. 

Finally, let hsr be the SR-only mobility (computed for 
strictly parallel interfacesji^ i.e. 5 = 0) and /.Jsr+rsr the 
total SR-f RSR mobility. Let us define the RSR mobility 
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FIG. 1. Surface roughness profiles computed for different g’s 
(A = 5 = 0.47 nm and tin = 2 nm). Note that the Si02/Hf02 
interface is actually smooth when A = 5 and g = — 1. The 
g — —0.5 stack approximately lies on the curve “A' = 0.47 
nm, A = 2.6 nm” of Fig. |4l 

A^rsr “ A*sr-i-rsr ~ Msr- Then, we expect from Eq. 
([S]) in a semi-classical Kubo-Greenwood approach: 

^RSR — “C'RSRA5g, (6) 

where Crsr is a constant. The RSR mobility is therefore 
inversely proportional to the rms fluctuations A and 5, 
and to the overlap g. Note that /trsr (as defined above) 
can be positive or negative, i.e. remote surface roughness 
can either enhance or inhibit the flow of carriers depend¬ 
ing on g. We will discuss the underlying physics later. 

In order to verify these assertions, we have performed 
Non-Equilibrium Green’s Functions (NEGF) calculations 
of the mobility in nFDSOI devices along the lines of Ref. 
[13 (same structures). We have considered tsi = 4 and 
tsi = 8 nm thick (001) Si films (dielectric constant esi = 
11.7) on a 25 nm thick buried oxide and a grounded n- 
doped Si substrate. The front gate stack is made of an IL 
of Si02 with varying thickness (eil = 3.9), a 2.4 nm 
thick layer of Hf02 (eHf 02 = 22), and a metal gate. We 
use the effective mass approximation for electrons. The 
wave functions can penetrate in the IL. Electron-phonon 
(PH) interactions are included in all calculations, but 
cancel out in /trsr (see Ref. [H for the methodology and 
a discussion about Matthiessen’s rule). 

The advantage of NEGF is that it makes use of an 
explicit, real space description of SR/RSR disorders and 
screening (no need for an approximation to the very com¬ 
plex Rr as in Kubo-Greenwood calculations). This al¬ 
lows for accurate calculations of SR mobilities.— In the 
present case, we have computed mobilities in samples of 
SR/RSR disorder satisfying Eqs. ([3]) (see Fig. [^, and 
averaged over a few configurations until convergence. 

As an illustration, the total SR-I-RSR mobility 
/ig^_i_RgR = -I- TJrsr is plotted as a function of carrier 

density n and g in Fig. [3^, and as a function of 5 and g 
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isi = 4 nm, ijL = 2 nm, A = 0.47 nm 



(b) n = 10^^ cm ^ 



FIG. 2. (a) /iSR+RSR as a function of carrier density n and 
g in a nFDSOI device with fsi = 4 nm and tm = 2 nm 
{A = 5 = 0.47 nm) (b) Msr+rsr (10“^ V.s/cm^) as a function 
of S (nm) and g in the same device, at n = 10^® cm“^. The 
red dots are the NEGF data, the blue surface is the best fit 

f^SR+RSR = f'SR + C'j^SR^'^5' 


in Fig. [5 }d. As expected from Eq. /^SR+RSR is pretty 
linear with 6 and g, and does not show, in particular, 
very significant oc contributions from the uncorrelated 
term in Eq. ©• We have also checked that RSR at the 
Hf 02 /Metal gate interface had much less impact on the 
mobility, at least for the Hf 02 thickness > 2 nm consid¬ 
ered here. 

The RSR mobility computed for g = —1 is plotted as a 
function of carrier density for different film and IL thick¬ 
nesses in Fig. |3l Since ^rsr > 0 when g < 0, the total 
mobility decreases when the fluctuations of the IL thick¬ 
ness are anti-correlated to the SR profile. Indeed, the 
carriers, squeezed by the surface electric field Eg, tend to 
localize where A(r) > 0; localization and scattering are 
enhanced by the concomitant increase of Eg when g < 0 
(IL tends to thin where A(r) > 0), and are, conversely, 
inhibited by the concomitant decrease of Eg when g > 0 



FIG. 3. RSR mobilities in 4 nm thick (dashed color lines with 
symbols) and 8 nm thick (solid color lines with symbols) Si 
films with A = 0.47 nm and g = —1. Black lines are Eqs. (0 

0 . 


(IL tends to thicken where A(r) >0). As expected, the 
RSR mobility decreases with increasing carrier density 
(or equivalently Eg) as the carriers move closer to the 
top gate stack; Hence RSR is most efficient (as is surface 
roughness) in the strong inversion regime, at variance 
with RCS for example. As such, it is weakly dependent 
on the thickness of the silicon film in the investigated 
range tsi > 4 nm. It shows, though, a strong (exponen¬ 
tial) dependence on the thickness of the IL; the thinner 
the IL the stronger the RSR correction. 

From these calculations, we can extract the prefactor 
C'rsr(^)Hl) as a function of the carrier density n and 
thickness of the IL. The Fourier components of the elec¬ 
tric field suggestjid^ that ^ Expanding 

K and a in powers of n, we tentatively fit the NEGF 
data with 


^RSR ~ 



exp 



tlL 

to 


( 7 ) 


and get a = 2.96 x 10“^ V.s/cm^, ni = 9.6 x 10^^ cm“^, 
to = 1.87 nm, and = 4.1 x 10^^ cm“^. Eq. ([7]) repro¬ 
duces the NEGE data very well for n > 2 x 10^^ cm“^, 
as shown on Eig. [S] 

How large can RSR scattering be in actual devices 
? As an illustration, we make the following assumption 
for the correlations between the Si/Si 02 and Si 02 /Hf 02 
interfaces 


/i(Hl) = , (8) 

where A is a vertical correlation length. The total 
PH+SR+RSR-bRGS mobility is plotted a function of tm 
in Eig. 131 for a 8 nm thick film in the strong inversion 
regime, and for different A and A'. RGS was modeled 
as a distribution of positive charges at the Si 02 /Hf 02 
interfacei^d^ with density urcs = 3.5 x 10^^ cm“^. In the 
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FIG. 4. Total mobility (PH+SR+RSR+RCS) in a 8 nm thick 
nFDSOI film as a function of tm, for different rms roughness 
A' at the Si02/Hf02 interface, and for different vertical cor¬ 
relation lengths A [Eq. (|8|] (A = 0.47 nm, urcs = 3.5 x 10^® 
cm~^). The carrier density is n = 10^® cm“^ (effective field 
Eeff ~ 0.8 MV/cm). The data were interpolated from a set of 
NEGF calculations. The lines are dotted wherever S > tm/S, 
i.e. when the IL might be pierced in many places. 

absence of RSR, the mobility is weakly dependent on the 
IL thickness, at variance with experimental datai^ The 
RSR can, however, significantly enhance scattering in 
thin ILs. It is, notably, very detrimental if the Si 02 /Hf 02 
interface is smoother than the Si/Si 02 interface, because 
the IL thickness fluctuations become anti-correlated to 
the Si/Si 02 roughness profile ((5 —^ A and 5 —>■ — 1 when 
A' —>■ 0). Note that alloying Si 02 with Hf 02 near the 
IL/High-K interface might blur the variations of the di¬ 
electric constant and effectively smooth that interface. A 
more detailed comparison with experiments^ would call 
for a comprehensive de-embedding of the different scat¬ 
tering mechanisms^iS Yet we can conclude that RSR, 
together with RPHSi^ (not accounted for at present in 
NEGF), are serious candidates in explaining the trends 
measured in very thin IL devices at strong inversion. 

Finally, we would like to discuss the effects of non¬ 
parabolic (NP) corrections and exchange-correlation 
(XC) effects. NP corrections can be included using a 
two bands k • p modeliSiSi for the conduction band, 
yet at a much larger numerical cost. Long-range XC 
effects have been modeled by an image charge self¬ 
energy correction (numerically computed along the lines 
of Ref. [H) , and short-range XC effects by a local density 
approximation ! ^^ 1 ^^ Calculations on a few representative 
configurations show that NP and XC corrections both 
decrease the mobility and Crsr (be., enhance RSR scat¬ 
tering), by a total of ^ 20 %. 

To conclude, we have discussed the impact of remote 
surface roughness scattering (RSR) at the Si 02 /Hf 02 in¬ 
terface on the electron mobility in FDSOI thin films. We 
show that the RSR corrections are dominated by the 


cross-correlations between the two interfaces of the IL. 
The mobility increases when the Si02 thickness fluctua¬ 
tions are correlated with those of the Si/Si 02 interface, 
and decreases otherwise. SR and RSR can not, therefore, 
be described as two independent mechanisms, and must 
be modeled concurrently. The RSR mobility shows an ex¬ 
ponential dependence on the IL thickness. In particular, 
RSR might be significant at strong inversion in thin IL 
devices where the Si 02 /Hf 02 interface is (dielectrically) 
smoother than the Si/Si 02 interface. 
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